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On the Feasibility of Real-Time 
Prediction of Aircraft Carrier Motion at Sea 

MENAHEM SIDAR^ 

BRIAN F. DOOLIN 

Ames Research Center, NASA 

Moffett Field, Calif. 9403S 

Abstract 

Landing aircraft on board carriers is a most delicate phase of flight 
operations at sea. The ability to predict the aircraft carrier's motion 
over an interval of several seconds within reasonable error bounds may 
allow an improvement in touchdown dispersion and a more certain value 
for ramp clearance due to a smoother aircraft trajectory. Also, 
improved information to the Landing Signal Officer should decrease the 
number of waveoffs substantially. 

This paper indicates and shows quantitatively that, based on the 
power density spectrum data for pitch and heave measured for various 
ships and sea conditions, the motion can be predicted well for up to 
15 seconds. Moreover, the zero crossover times for both pitch and 
heave motions can be predicted with impressive accuracy. 

The predictor was designed on the basis of Kalman's optimum 
filtering theory for the discrete time case, adapted for real-time 
digital computer operation. 

^NRC Research Associate 


1. Introduction 


The landing phase of an aircraft aboard an aircraft carrier 
represents a complex operation and a demanding task. The last 10 to 15 
seconds before aircraft touchdown involves terminal guidance and 
control problems, where not only the aircraft is disturbed by several 
kinds of stochastic (wind) disturbances, but also the touchdown point 
(on the ship) is being moved randomly. Despite the wind disturbances 
and the final point (target) random motion, the landing accuracy 
specified for carrier operations is very high, i.e., a few tens of 
feet longitudinal landing dispersion. Such a terminal point problem 
is made tractable in a most natural way by assuming that the ship's 
position can be predicted for several seconds ahead so that the airplane 
is guided toward the future position of the touchdown point. The scope 
of this study was to establish to what extent a stochastic process, 
like the ship's motion, is predictable over moderate periods of time. 

Quantitative results obtained throughout this predictor's 
feasibility study concerning the relationship between the prediction 
error versus the prediction time, the influence of measurement noise, 
and the "narrowness" effect of the ship motion power density spectrum 
are presented. Digital simulations show that the prediction accuracy 
does not degrade prohibitively even for quite large measurement noise. 

A variety of possibilities with respect to the incorporation of 
the prediction algorithm in the Aircraft Carrier Landing System (ACLS) 



can be investigated, but those topics are out of the scope of this 
paper. Their common denominator consists of the ability to predict 
the ship's motion in heave and pitch for periods of 12 to 15 seconds. 
Feasibility of predicting the ship's motion within acceptable bounds 
of error can also lead to improvement of the Landing Signal Officer 
(LSO) decision policy for waveoffs. 

The need for prediction for carrier landing operations was 
pointed out several years ago by Durand [1], Durand and Wasicko [2], 
Kaplan [3], and Siewert and A'Harrah [4]. Loeb [5] also indicated 
the need for predicting the ship's motion. 

A tremendous amount of theoretical modeling, experimental results, 
and the collection of a large amount of data over the years are 
discussed by Powell and Theoclitus [6], Kaplan [3], Johnson [7], and 
many others. 

A study of prediction techniques for aircraft carrier motions at 
sea was done by Kaplan [3] , who considered a deterministic technique 
based on a convolution integral representation with wave height measure- 
ments at the bow serving as input. He derived, from the ship response 
time-history functions, a Kernel -type weighting function which operated 
on the measurements in order to provide the predicted motion history. 
Model test data indicated that this technique yielded ship's pitch 
prediction for up to 6 seconds, but the method suffers from severe 
limitations and practical implementation difficulties. A hybrid pre- 
diction technique, based on modern control theory, was suggested in [jj 
as a possible future approach. 
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Our approach was to make a rather direct use of the ship's motion 
characteristics and the measuring instrumentation existing on board the 
ship in order to get the predicted motion. Using this information, a 
predictor based on Kalman's theory of optimum estimation was designed. 

Several circumstances contribute to the success of this approach. 

The size and mass of the ship significantly filter the motion of the sea. 

A complete landing operation is short enough that the stochastic processes 
are reasonably taken to be stationary. Finally, the prediction interval 
is only a small fraction of the time it takes each aircraft to land. 

This paper is divided into three parts: 1) the derivation of the 

mathematical model of the ship's motion, 2) the rationale for the 
predictor implementation, including the Kalman filter and predictor 
equations, and 3) discussion of some of the results obtained. Since we 
are interested only in the most critical aspects of the landing 
operation, namely the characteristics of the longitudinal channel, we 
merely investigate in the sequel the predictability of the pitch and 
heave motion of the ship. 

II. Modeling the Carrier Motion 

As we pointed out before, quite a large amount of data exists 
describing the motion of aircraft carriers at sea. Extensive experi- 
ments and simulations have been carried out, both at sea and in water 
tanks, establishing frequency response curves and power spectral 
density functions (psdf) as a means of representing ship pitch and 
heave motion characteristics for utilization in a systems analysis of 
the whole carrier-aircraft landing system. 
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Fig 1 
Fig 2 


Power spectral density functions describing, globally, the 
statistical behavior of the pitch and the heave motions have been 
established and measured for several types of carriers at different 
sea conditions [6], [7], An analysis and a close comparison of these 
data reveal that the psdf, <t(w), where w is the carrier motion frequency 
in radians per second, is not affected too sensibly either by the type 
of the carrier or by the sea conditions (see [7]). Moreover, the 
function peaks sharply around a center frequency of about 

= 0.60 rad sec"^ (Fig* 1)» the sea state changing the value of the 
peak at u) = 0 )^ (see [7]). 

Data obtained from basin model experimentation confirm (Fig. 2) 
the measured power density spectrum, showing the same narrow-band 
aspect, but centered at = 0.75 rad sec’^. 

It is obvious, also, that for the prediction periods of interest, 
we can make the plausible assumption that the ship heave motion £(^) 
and pitch motion 6(^) are stationary, narrow-band, stochastic processes. 
Both processes are actually continuously measured aboard the ship as a 
part of the SPN-42 ACLS, the measurements being contaminated by random 
noises v(£) and w (t^) , respectively (see [8]). 

In order to obtain the mathematical model of ♦^(to) and ‘+’g(t*i) for 
the ship's heave and pitch, we take the innovation process point of 
view; namely, we assume that z(t) and 9(t) are stochastic processes 
generated by a white, Gaussian, random process passing through a causal 
and invertible lumped transfer function G^(s) or Gg(s) (see [9]). 

The first step of the procedure is to approximate the experi- 
mentally obtained density (unctions and (<n) by analytic 
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i expressions that accurately represent the important part of the spectrum. 

I 

K 

t As a result, one obtains a single, equivalent transfer function 

i * 

I G (s) = G„(s) = G (s) of the filter acting on the Gaussian noise, which 

I — z — O “ — S ” 

I is essentially the same for both the pitch and heave motions and is 

i 

I 

f given in (1): 

* I! 

I “ 

? 

with the following nominal values for g, a, and d: g = 0.6, a = 0.06, 

d = 0.36. This approximation is "pessimistic" with respect to the 
prediction ;^>roblem, as is shown by the comparison with 4>(u») in Fig. 1. 

The reason for this pessimistic approach is twofold; 1) to avoid running 
into high-dimensional systems unnecessarily, and 2) to obtain conserva- 
{ tive values for the maximum achievable prediction time. 

From (1), by making use of the usual "half power" definition of 
the bandwidth for a narrow-band process, equivalent bandwidth (BW) of 
the ship's motion process transfer function G^(id) is obtained: 

BW = ^/l + 2C - /I - rad sec'^ centered at (2) 

For = 0.6 rad sec"^ and C = 0 . 05 , BW = 0.1 * 0,06 rad sec"^. 

Idealized narrow-band random processes, also called ideal bandpass 
stochastic processes, have been studied by Rice [10], who obtained 
valuable theoretical results. In particular, the autocorrelation 
formula and the probability distrilxition of zero crossing for the ideal 
bandpass process are given. In this paper we will validate those results, 
comparing them for our particular problem of ship motion predictability. 

(i- 



Once 9s obtained, It is necessary to express it 
equivalently in the state space form. The equations take the form of 
a pair of linear differential equations driven by a Gaussian white 
noise process. The s in the numerator of (1) normally implies that 
the derivative of the input function u(t) is to be used as a forcing 
term in the state space representation. To avoid differentiating the 
random process, an adequate transformation was performed [11] with the 
following set of differential equations being obtained, where x(t) repre- 
sents either pitch or heave, u(t) is the scalar random input, v(t) is 
scalar noise in the measurements, and the dot indicates the time derivative 


? * ^ 

(3a) 

y * c-x + V. 

(3b) 

In (3), the following is obtained from (1): 


A r 0 n 

^ • L-^ -J 

(4a) 

b- = [g, -ga] 

(4b) 

J « [1. 01. 

(4c) 

The following assumptions are made with respect 

to the noise and 


the input: 



E [u] = E [v] = E 

[v x.] = 0 

(5a) 

! ' ? 

• 6{t-s) 

(5b) 

E (v(t) v(5)l = R 

• 6(t-s). 

(5c) 
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Fig 3 


The system of equations (3) is controllable and observable. In these 

T 

equations, x^(t) in x> > (x^, X 2 ) represents position, either heave or 
pitch angle. Equations (Sa) through (5c) express the assumption that 
the ship motion and the measurement noise are uncorrelated. 

III. The Predictor Equations 

Several approaches and algorithms are available for obtaining 
x(t+T), the prediction of x(t) at t>0 seconds frwn now: 

1) the Wiener approach (12), [13], which assumes stationarity, a 
plausible assumption over the short periods of time required for landing; 

2} the Ragazzini-Zadeh approach [13], [14], which is useful for 
the finite>time measurement case, but is otherwise complex; 

3) the time series analysis and prediction algorithm [IS]; 

4) the Kalman predictor approach [16], [17]. 

The last approach is adopted here because, in avoiding cumbersome 
computations, it makes real-time digital computation possible. 

The best estimate of the ship's motion at time t, te [o*^f] 
is denoted by x(t). It is the conditional expectation of x(t) based on 
all prior measurements y(t). Then the linear least-squares prediction 
theory of Kalman [16], [17] gives for the best predicted motion, 
x(t+T), the expression: 

X (t + T) * <j) (t + T, t) X(t) (6) 

where 4(t,o) is the transition matrix for (3). 

The ccHsputation process for x(t ♦ t), therefore, divides into the 
following two steps: 1) calculate x(t); then 2) use (6) to obtain 

x(t ♦ t) for the desired prediction t (see Fig. 3). 
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Since the Kalman formulation was adopted to enable real-time 
digital computation, (3) - (6) should be replaced by their discrete 
form. The discrete-time representation of the ship motion is given by: 


^ • s- ik * ik 

where: 


(7a) 

(7b) 



> a sequence of random uncorrelated measurement noise: 

■ a sequence of random uncorrelated inputs: E(uj^u^) ■ 

■ the sampling time. 






The discrete-time version of the optimum filter is obtained from the 
following set of equations: 


ik ■ *k,k-i 5k- 1 * !Jk ♦k.k-i ?k-i> 

h ■ h •!;■ -is' V * • 


(8a) 

(8b) 


Note that the measurement vector ”c” has been replaced for convenience 
by H matrix "C” where: Cj^ » (J 

The a posteriori covariance matrix is obtained from: 


-k * -k-1 * -k-1 


(8c) 


- 9 - 



and the a priori covariance matrix is given by: 


?k ■ fk' ■ !Jk ? L’k'- 


{Sd) 


-k* -k ' -k Kalman filter gain) are square matrices. By 
time k, as a result of (8a) through (8d), the best estimate Xj^; 




’’past" values of yj^I. h e (0,N) is generated, 


The initial value of the error convariance matrix Pj^ is given by 


the matrix: 


Eo * io !o- 


(8e) 


which is used as a startup value for the recursive scheme. 

In this specific case we have to compute only two optimal gains, 
namely K , and K,, , given by: 




?iik> 
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Bilk’- 


(9) 


The predicted vector x(t*T) is obtain^ fr«a (10) which is the 
discrete equivalent of (6): 

?k*B " * ik 


( 10 ) 


where m » ^/Ts‘ 

The transfer matrix i^(k*m,k) is given by: 
(k+m,k) ■ ^ (t) « 

Fcos 8t ♦ (a/26) sin 6t 
I (-d/B) sin 6t 


(1/6) sin Bt 
cos Bi - (a/ 26 ) sin 



( 11 ) 


and 6^ * d - a^/4. 
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Fig 4 I 5 


Finally, since we are interested only in x^Ct^t), the ship's 
optiaal predicted value for t seconds ahead, one obtains the following 
result: 

5 |(t+T) » (cos St (a/2S) sin 5t) Xj (t) 

+ (1/8) sin 6t ‘X^ (t)]. (12) 

For a narrow-band process, which is the present case, the parameter a 
in (12) is very small. Thus, for moderate prediction times, (12) can 
be ipproximated by: 

Xj(t+T) * Xj(t) cos St ♦ x^lt) (sin 6t)/6 . (13) 

Equation (13), which will be used later in calculating the auto- 
correlation function, shows that under the stated conditions, the extra- 
polation is equivalent to predicting the sta* of an harmonic oscillatur 
on the basis of an estimate of its present state. That this approx i mat !■ 
is appropriate is easy to see from the quality of the prediction, resuits 
of which will be discussed next. 

IV. Results and Discussion 

The data shown in Figs. 4 and 5 sinnaarize the results of a large 
number of cases investigated by a digital simulation of the ship’s mutiur 
aiwl its prediction. The figures show the effect of the prediction time 
T in terns of two performance criteria of significance in the carrier 
landing problem. 

The first, J^, measures performance in terms of the error in 
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predicting position. The error in prediction, e^, is defined as; 
Ep (t -n) = Xj(t +t) - Xj(t +t) , “V't c [0,:] 


or 


(kT^+T) ixj(kT^+T) - Xj(kT^+T) , -V' k e [0,N], (14) 


The criterion is defined as the following scalar functional: 


1 ^ r -v-1 

Jl = (a ) 


IT 


/a/N+1) I e: (kT +T). 
k=0 C "! 


(15) 


The scalar depends on the spectral width a and the measurement noise 
covariance R, in addition to the prediction time t, as Fig. 4 shows. 
There it can be seen that if the measurement noise is not prohibitive, 
prediction times as long as 10 to 15 seconds are attainable with 

2 

reasonable accuracy. For example, a 5% measurement noise (R = 0.02o ), 

^1 

with a =1 meter (3.1 feet) and t = 10 seconds, gives 
-1 

a - a J. z 0.24 meters (0.8 feet), which can be considered acceptable. 

PI 

The results obtained during the study suggested the introduction 
of a second criterion with respect to the predictor quality. It 
covers another important practical use for the predictor, the prediction 
of crossover times. This second criterion is defined as follows: 


J, ^ (o) /2n) 
-2 0 


M 

I 

i=l 


^ -CRO 




(16) 
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where M is the number of crossover points counted in the fixed, finite 
interval of an experiment, and 


‘ICRO ‘ (Ti:x(t.»,).01 


(Tj:x(t^+At^+T)*Ol 


( 17 ) 


where 


0 i «! < - t. 


Here T^:x(tj^+t) means the ith time that x crosses zero. Thus is 

defined as the difference between the crossover time of Xj(t+x) and the 


nearest successive crossover time of Xj(t+t). It tells how long the 
actual ship heave, for example, differs in sign from that predicted. 

Fig. 5 shows that this error is small, even for prediction times of 
10 to 15 seconds. For example, under the same conditions just given, 
namely a 5% noise-to-signal ratio and a 10-second prediction time, and 
for seconds, the total time that the predicted motion 

differs in sign from the actual motion is 0.5 seconds over a measurement 
interval of 60 seconds. The practical implication of the results shown 
in Figs. 4 and 5 is clear and important: both the motion of aircraft 
carriers at sea and the crossover times of their motion can be predicted 
accurately over periods of time that are long enough to be operationally 
useful. 

These results are in good agreement with Rice's analysis [10] of 
ideal narrow bandpass processes. However, the comparison has to be 
made with care since in our study the noise enters into the measurement? 
as well as at the input. Furthermore, our results are derived from 
simulations of finite duration. Nevertheless, the fact that our 
results are consistent with his analysis of narrow-band noise 
strengthens our confidence in the predictability of the motion. 
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Evidence that the predictability of the motion of the carrier is 
due to its effective narrow-band character is seen from comparing the 
high correlation of our results with that to be expected of a narrow-band 
process. From eqn. 3.2-5 in [10] one gets the following expression for 
the autocorrelation function R(t) of such a process: 

R(t) = o^ [sin (■« BW) T / t] cos t (18) 

Setting T equal to an integral number of periods, where f^=l/T^ 

is the center frequency of the bandpass (or the peak frequency of the 
power spectrum of the carrier's motion), one obtains: 

R(t)/o^ * 0.2 IT ^ / 0.2 TT (19) 

This gives R(T^) = 0.935a^, and R(2T^) = 0.757o^. 

These results compare well with ours; the correlation between the 
predicted values of x^(t), namely (t+x), and their actual values is 
very nearly unity over times of the order of 10 to 15 seconds. 

Furthermore, the influence of the process narrowness BW can also 
be obtained from (18) and compared with the results in Fig. 4. The 
autocorrelation function R(x) is a sine function (see 18) of the 
bandwidth, a result which is in agreement with the results obtained in 
our study. 

Rice also made an analysis of the expected zero crossings for a 
narrow-band process, and again, interpreting his formulas in our terms 
gives values consonant with the results shown in Fig. 5. 

Eqn. 3.3-12 of reference [10] gives the expected number of zero 
crossings of Xj(t) per second, as: 

. 2 [ 1/3 - fj’) / (f^ - f^) (20) 
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;'!|iP!f!!»;i'iK!»Wilii:i|ip!lt!ll||iKW^^^ 









where and are the upper and lower frequency limits for the ideal, 
band limited process, and BW * ®PP^®*ches f^, as in 

the narrow- band case. 


N 


-2 




= 2 f 


-o 


( 21 ) 


which is the number of zero crossings for a sinusoidal heave (or pitch) 
motion with frequency f^=u)^/2ir. This result agrees with ours (Fig. 5), 
even for prediction times as large as twice the period T^. 

Also from [10], the probability that a second zero crossing of 
Xj(t) lies within t and t+At, at a particular time t, say 
t a K6 = K (Tq/ 2), is given (approximately) by (22) for the ideal 
narrow-band process: 


pr (n ) = (1/2) (Y / [1 + (t-ke)^]^^^) 

-At 


( 22 ) 

t = k(T^/2) + At 
k a 1,2,... 


where 


Y = /3 [(fj^ + fj^)^ / BW] ; 6^1/ (f^^ + f^ . 

Equation (22) shows that the probability density is a symmetrical and 
decreasing function around t = k(T^/2), its peak value and dispersion 
depending on the bandwidth BW. For example, in the specific case of 
the ship's motion prediction, with BW = 0.1 f^ and = 10.45 seconds, 
pr(n ) = 3.46, whereas pr(n ) * 0.01. This means that the 

■' ‘5o 'h 

probability of a zero crossing at (T^ +1.0 second) is very small 
(0.01). Here again, our results show a nice regularity of the 
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zero-crossing points which are (almost) identifiable with the 
zero-crossing points of an harmonic ship^s motion with frequency f^. 
This is, basically, the reason why it is possible to obtain a good 
prediction of the zero-crossing events for Xj(t+x), even for relatively 
high values of t. 

V. Conclusions 

The feasibility of predicting aircraft carrier motion at sea by 
measuring the actual ship's position was investigated. The ship's 
motion mathematical model based on statistical data, such as power 
density spectrum representation, was established. Subsequently, a 
discrete-time Kalman filter-predictor adapted for real-time computation 
on a digital computer, was investigated. The results obtained show 
that a maximum achievable prediction time of up to 15 seconds can be 
reached within reasonable acceptable errors. 

Being able to predict accurately the ship's motion can lead to an 
improvement of aircraft landing accuracy. This can be accomplished, 
for instance, by generating new terminal guidance (landing) laws 
making use of the future ship's position. Moreover, the possibility 
of prediction can eventually improve the LSO information and policy 
for landing acceptance or wave-off. The possibility of processing the 
measured motion by Fast Fourier Transform Algorithms (FFT), in order to 
obtain the power density spectrum for the ship's motion in real time, 
may lead toward an adaptive predictor scheme. 
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Fig. 1. Ship motion heave spectrum. 
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Fig. 2. Pitch and heave representative spectrum. 
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